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ABSTRACT 

Differential-phase partial-reflection measurements have been used to deduce 
D-region electron -density profiles. The phase difference was measured by taking 
sums and differences of amplitudes received on an array of crossed dipoles. The 
reflection model used was derived from Fresnel reflection theory. Seven profiles 
obtained over the period from October 13, 1971 to November 5, 1971, are presented 
along with the results from simultaneous measurements of differential absorption. 
Some possible sources of error as well as error propagation are discussed. A 
collision-frequency profile is deduced from the electron concentration calculated 
from differential phase and differential absorption. It is concluded that the 
differential phase measurement would be a useful addition to the partial -reflec- 
tion experiment. 
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1. INTRODUCTION 

The D region (50-90 kilometers) is probably the most complex and least under- 
stood region in the ionosphere. Efforts to describe the D-region have been ham- 
pered by the complexity of the neutral and ion chemistry and the difficulties 
associated with measurements at these altitides. The parameters of the lower 
ionosphere which can be used to understand the processes that occur must be mea- 
sured accurately and reliably. One parameter, which is important in the iono- 
sphere and particularly in the D region, is the electron concentration. One 
goal of an experimental aeronomer, then, is to develop a system which can measure 
the electron concentration 'with sufficient accuracy that a reasonable description 
of the ionosphere can be made. 

Measurements of electron concentration can be made by sounding rockets or 
by ground-based techniques. Rocket measurements of D-region electron concen- 
trations and collision frequencies have become quite sophisticated and are be- 
lieved to be accurate (Mechtly, e£ al_. 1967). However, because of cost limita- 
tions, only a limited number of rockets can be used to gather information. 

Also, sounding rockets can only be launched from relatively few locations. The 
disadvantage of ground-based experiments is that the data are more difficult 
to interpret. However, ground-based experiments can be set up in a variety of 
locations and can be operated continuously or as required. Thus, large amounts 
of data from several locations can be used to study diurnal, seasonal, and geo- 
graphical variations. 

One such ground-based technique for measuring electron concentrations in 
the D region is the partial -reflection experiment. Weak reflections from the 
D region are used to calculate the electron concentration over a range of alti- 
tudes. Gardner and Pawsey (1953) were the first to use a partial -reflection 
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system to obtain electron -concentration profiles. Since then several experi- 
menters, for example, Belrose and Burke (1964), Pimat and Bowhill (1966), 

Thrane, et al . (1967), Gregory and Manson (1967), and Birley and Sechrist (1971), 
have developed the partial -reflection experiment to the point where it provides 
useful quantitative measurements of the electron -density distribution in the 
lower ionosphere. 

Most of the partial -reflection measurements that have been made are of 
differential absorption. Fejer (1961) suggested that the partial -reflection tech- 
nique be modified to include measurements of the phase difference between the two 
magneto-ionic modes. Differential phase measurements would add a simultaneous 
independent electron concentration measurement which along with differential 
absorption might subsequently be used to deduce a collision -frequency profile. 
Belrose (1970) also suggested the measurement of differential phase to improve 
the partial -reflection technique. 

Two methods of measuring the phase of a received signal are possible. The 
first method, suggested separately by Connolly (1971) and by Austin (1971), com- 
pares the phase of each magneto-ionic mode in order to obtain a voltage propor- 
tional to the difference in phases. Connolly did some calculations of expected 
experimental results. He concluded that differential phase measurements when 
combined with differential absorption would provide better electron density and 
collision frequency data from 70 to 90 kilometers. Austin performed a direct 
phase measurement and obtained one electron density which showed reasonable val- 
ues for electron concentration. The experimental arrangement that he used is 
described by Austin, et al . (1969). 

The second method which was used by von Biel, et^ al_. (1970) and von Biel 
(1971) obtains the phase difference by measuring mean-squared amplitudes with 
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different antenna polarizations. From these amplitudes a correlation coeffi- 
cient is calculated which is dependent on the integral over the propagation 
path of the phases of the two magneto-ionic modes. These investigators obtained 
reasonable results, although extensive data collection was not made. 

At the University of Illinois Aeronomy Field Station a system for measuring 
differential phase has been set up in conjunction with the existing partial -re- 
flection experiment. The method used is similar to that of von Biel, et al . 
(1970) and von Biel (1971). In Chapter 2 definitions and relations necessary to 
calculate electron concentrations from differential phase measurements are given. 
A development of the Fresnel reflection theory along with volume scattering the- 
ory is presented in Chapter 3. Also a discussion of the significance of each 
theory as applied to differential phase partial -reflection measurements is given. 
Chapter 4 describes the system used to obtain phase measurements as well as the 
procedure for making these measurements. Electron concentration profiles ob- 
tained by this system are presented in Chapter 5. A critical evaluation of the 
system is made in Chapter 6. Possible sources of error and their effects on 
electron concentration calculations are considered. Also, the differential 
phase system is evaluated as a possible tool for studying D-region ionization. 

In the final chapter conclusions and suggestions for further work are given. 

The computer programs used to process the differential phase data are listed in 
the appendix. 
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2. RADIO WAVES IN THE IONOSPHERE 

An electromagnetic wave propagating through a medium can be described by 
the expression for the electric field intensity vector 


? = E exp[j(wt - lc • r)]. (2.1) 

where tu is the operating angular frequency, t. is the propagation vector in the 
direction of phase propagation, and r is the position vector. As the wave pro- 
pagates through the medium, the magnitude and phase of ti change as lc and r 
change. These changes in lc imply that the medium has an effect on the wave as 
it propagates. By measuring the changes in the magnitude and phase of infor- 
mation about the medium such as the electron concentration can be deduced. 

The medium which will be considered here is the D region of the ionosphere. 
Since the earth's ionosphere consists of a homogeneous admixture of neutral and 
ionized species immersed in the earth's magnetic field, it is termed a magneto- 
ionic medium. As a magneto-ionic medium, however, the ionosphere is relatively 

weakly ionized. The concentration of free electrons which are largely produced 

7 9-3 

by photoionization is on the order of 10 to 10 m during the day. Both the 
magnetic field and the ionization have a significant effect on the propagation 
vector k. 

Another important parameter of a magneto-ionic medium is the collision fre- 
quency. The number of collisions between electrons and neutral particles per 
second is termed the collision frequency and is represented by the symbol v^. 
Other collisions such as electron-electron, neutral -neutral , or neutral-ion do 
not have a significant effect on radio waves due to the mass of the particles 
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or the small number of collisions. A knowledge of the electron-neutral collision 
frequency as a function of altitude is necessary for the expressions which des- 
cribe the effect of the medium on the wave. A more detailed treatment of colli- 
sion frequencies and their measurement is given by Lodato and Mechtly (1971). 

The propagation vector £ depends on the refractive index n of the medium 
and is given by 



( 2 . 2 ) 


where X is the wavelength of the operating frequency and k is a unit vector. 

The refractive index is defined as the ratio of the speed of light in a vacuum 
c to the phase velocity of the wave v^ in the m gneto-ionic medium and is given 
by 



(2.3) 


where is the relative permittivity of the medium. The relative permeability 
will be taken to be unity for the ionosphere. 

A wave traveling in the z direction will have components of the electric 
field intensity vector and the magnetic field intensity vector H in the x and 
y directions. The polarization of the radio wave is defined by 


R = E /E 
y x 


= - H /H 
x y 


(2.4) 


When the wave polarization R is equal to a constant the wave is propagated 
through the medium without changing its polarization. For a magneto-ionic 
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medium there are two of these characteristic waves, the ordinary and the extra- 
ordinary. The absorption and phase shift of the extraordinary wave is signifi- 
cantly greater than for the ordinary wave. North of the magnetic equator the 
field vectors of the ordinary wave rotate from north into west, and the field 
vectors of the extraordinary wave rotate from north into east. 

The effect of the magneto-ionic medium on R can be described by consider- 
ing a wave traveling through the medium. 


t = £ exp [ - j (lc • r) ] 


(2.5) 


The time dependence has been suppressed. The propagation vector in one direc- 
tion (z) is given by 



(2ir/X) 


| n dz 

h °’ X 


( 2 . 6 ) 


The subscripts o and x refer to the ordinary and extraordinary waves, respectively. 
By substituting Equation (2.6) into (2.5), the expression for the wave can be 
written as 


E 

o,x 


£ ex p[-j (f— 


n dz) 1 

OX 

£ 


(2.7) 


By separating the refractive index into its real and imaginary parts (2.7) can be 
written as 


E 

o,x 


£ exp[- 


2it n* dz] • exp[-j(^- I n x dz)] 
o,x J r L J [ . o,x J J 


( 2 . 8 ) 
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where the superscripts * and i denote the real and imaginary parts of the re- 
fractive index. As the wave travels through the medium it is absorbed and the 
phase is changed. The absorption depends on the integral of the imaginary part 
of the refractive index over the path. The change in phase is a function of 
the integral of the real part of the refractive index over the path that the 
wave travels. 

The ratio of the returned electric field strength vector of the extraordi- 
nary wave to that of the ordinary wave is given by 


V2 0 = ex p[- 


2t r 


(n^ - n*)dz] • exp[-j 


. 2n 


( n x - n o) dz l 


(2.9) 


The phase angle of this ratio can be taken from equation 2.9 

(n r - n r )dz . (2.10) 

v o x ' 

Alternatively the expression for A<j> can be described as the difference between 
the phase angles of the ordinary and extraordinary mode of polarization. This 
difference A<f> is the phase difference which is to be measured. 

In the development of magneto-ionic theory an expression is derived for 
the refractive index in terms of the electron concentration, the collision fre- 
quency, the operating frequency, and the magnetic field. The classical Appleton- 
Hartree formulas were used by Gardner and Pawsey (1953) and other earlier experi- 
menters. A full treatment of the Appleton -Hart ree formulation is presented by 
Ratcliffe (1959). Sen and Wyller (1960) derived a generalized magneto-ionic 
theory and corresponding formula for the refractive index. The use of the 
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generalized formula is necessary for calculations in the D region. For mid- 
latitudes the quasi -longitudinal approximation can be used with negligible error. 
The quasi-longitudinal approximation to the complex refractive index is given by 


n = n 


% 

'V 


- 

3 n 

o,x 

J o,x 

2 


“N 

’0) ± 0)^ 

0)V 

V 

m 

L m 


JS 


3/2 


ru ± Co. 


m 


- j 5/2 


U)V 


m 


£ 


S/2 


r w ± 


m 


1/2 


( 2 . 11 ) 


where is the longitudinal component of the electron gyrofrequency. is the 
plasma frequency and is given by 


2 _ Ne^ 

me 

o 



( 2 . 12 ) 


where N is the electron concentration m ; e is the electronic charge C; m is the 
mass of an electron kg; and e Q is the dielectric constant Fm . p i s t ^e 

semiconductor integral. Approximate formulas for evaluating p are found in 
Hara (1963). The real and imaginary parts of the refractive index (n Q x and 
n* x ) can be found by expanding (2.11) by the binomial theorem. 



-j 1/2 


5/2 


(i)V 


m 


Ps/2 


0)1 u) Ll n 


m 



U)V 


m 




± 0 ), 


v J 
m ; 


- 1/2 


+••• (higher order terms) 


(2.13) 
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Since 


£ 

m 



3 

+1 . 

3 

, -V 

✓ ✓ 1 

r 1 “l 1 r 

3 

+1 

3 

5/2 

[ v 

" m 

X 

03V 

m ... 

v m ^ 3/2 
1 m J 

L m 1 


'v i 
% ] 


the higher order terms can be neglected, and the real and imaginary parts 
be equated. 


o,x 


1 - 


N 


03V 


m 


03 i 03. 1 


m 


3/2 


(JL) Z 0J T 'V 


m 


1 1/2 


J J 


2 

n 1 = 5/4 — £ - 

o,x 03V /° 5/2 [ v J 


( 2 , 


The real part can again be expanded by the binomial theorem. 

2 


03 . . 

r , N 

r + “ L i e 

'j 

3 

+ 

3 

Jl J. « 

o 2 oiv 

m 

1 V. j > 3/2 

1 V J 

L m ' 


( 2 . 


(l).. 

r . N 
n v = 1 - ~ — 

X 2 03V 


m 


r ■ “li >° 

•03 - 03^' 

v m 3/2 

v m ' 

• v m • 


The imaginary part of each mode is given by 

2 


n 1 = 5/4 


N >o 


03V 


m 


/° 5/2 


03 + 03. 


m 


( 2 . 


( 2 . 


n* = 5/4 — — P 

X 03V ' 

m 


“xi >=> 


5/2 l v m j 


can 


14) 


15 ) 


16 ) 


17) 


( 2 . 18 ) 
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If the change in phase can be measured over the path that the wave travels, 
the electron concentration can be calculated for the medium. For partial re- 
flection, however, the reflection process must be considered. Expressions for 
the electron concentration when reflection is included will be developed next. 
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3. THE THEORY OF PARTIAL REFLECTIONS 


The basis for the partial -reflection experiment is the reflection process 
in the D region. The incident wave is partially reflected as opposed to totally 
reflected. An understanding of the effects of the ionosphere on the wave as 
it passes through the ionosphere and as it is reflected is essential to the 
method . 

The nature of the reflections has been a subject of considerable study. 
Gardner and Pawsey used the Fresnel reflection theory which was also used by 
Belrose and Burke. The volume scattering theory of Booker (1959 ) was proposed 
by Flood (1968) as a better description of the reflection mechanism. Several 
investigators have studied the characteristics of partial reflections, for 
instance, Gregory and Vincent (1970) and Fraser and Vincent (1970). The 
heights of reflection and phase path changes, respectively, were investigated 
to gain further insight into the nature of the partial reflections. The two 
principal theories will be discussed here as they apply to differential phase 
measurements. 

3. 1 Fresnel Reflection Theory 

Assume a horizontally-stratified ionosphere where layers of different re- 
fractive indices constitute the stratifications. Furthermore, the refractive 
index changes from n to n + An in a distance small compared to the wavelength 
of operation. Therefore, the reflection coefficient from a boundary between 
two regions of refractive indices n^ and ^ is given by (Belrose and Burke, 1964) 



6n 

2n 


(3.1) 
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if 


2 



0) 


By differentiating the quasi-longitudinal approximation to the complex refractive 
index (2.11) with respect to N, the reflection coefficients for the ordinary and 
extraordinary waves are given by 


R = - 
x 


SNe 2 

*f“ + V 

1 7° 

3 

- + 

3 

2n me wv 

V I 

1 > 3/2 

V 1 

o o m L 

A m J 


1 m 1 

6Ne 2 

"i-J 

3 

i 

3 

1 & 

Zi 

3 

3 

2n me cov 

x o m - 

1 V j 

A m J 

I /° 3/2 

V J 

1 m 1 


j 5/2 £ 


+ 3 5 /2 


ru + a). 


5/2 


m 


CO - 0). 


(3.2) 


5/2 [ v m Jj 


The refractive indices are both close to one; therefore, n^ will be approximately 
equal to n Q . 

The part of the reflection coefficient which enters into the differential 
phase relations is the phase change due to reflection. The ordinary reflection 
coefficient can be expressed as 



(3.3) 


where y q is the phase shift that occurs due to reflection. By using equations 
(3.2) and (3.3), the phase y q can be written as 



(3.4) 
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Similarly, for the extraordinary wave the reflection coefficient is 


R 


x 




(3.5) 


where 


tan 



(3.6) 


The expressions for y q and are dependent on collision frequency and not 
on the electron concentration. 

For an electromagnetic wave that is transmitted vertically from the ground, 
reflected from a discontinuity, and received back at the ground. Equation (2.5) 
becomes 


^ = R ^ exp[-jlc • r] 


(3.7) 


where R is the reflection coefficient. Since the wave traverses the same path 
twice, the magnitude of the propagation vector in the £ direction becomes 



4trf 

c 


ndz . 

£ 


(3.8) 


By inserting the phase of the reflection coefficient for each magneto-ionic mode 
into (2.10), the total phase difference is given by 

_h 

( n o " n x )dZ + ^ Y x " To) • ( 3 - 9 ) 

0 
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If A<f> is measured at two different altitudes and h^, the difference be- 
tween the two measurements is given by 


A4>, - A4> = 


2oj 

c 


(n r - n r )dz + (y - y ) 
o x' v 'x 'o'2 


• Vl 


(3.10) 


where h^ > h^ and the subscripts refer to quantities at h^ and h^, respectively. 
Since the medium can be assumed to be nearly homogeneous between h^ and h ^ if 
the difference between h^ and h^ is small, the integral in (3.10) can be approxi- 
mated by 


4mf 

c 


rh. 


, r , - 4irf . r 

An dz % An 


Az 


(3.11) 


Az is the difference in heights, h 2 - h^ From the substitution of the differ- 
ence in the real parts of the refractive index given by 


An 


r 


r 

n 

o 



1/2 


Ne 

me oiv 

o 


ru) + (i). 


m 



(3.12) 


into Equation (3.10), an expression for the electron concentration as a function 
of collision frequency and differential phase can be written. The electron con- 
centration N is given by 


N = 


mce v 
o m 

2 

e Az 


[(A* 2 - A*p - (y x - y 0 ) 2 * (Y x - r Q J 1 ] 


~(b) - 0). ^ 

L y=> 

3 

« 

3 

V 

+ “Ll ^ 1 

'CO + 

tl «, J > 3/2 

1 

E 

v ) h 3/2 
m J ' 

V 

- m 


(3.13) 
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Expected values of the phase difference (Acf^-AtJ^) can be calculated using 
a model ionosphere. The electron-concentration profiles for rockets 14.244, 
14.246, and 14.247 (Figure 3.1) were used to calculate the reflection coeffi- 
cient difference Y 0 - Y x and the change in phase difference A<J» ^ -A<f> ^ . The corre- 
sponding collision -frequency profiles used for these calculations were given by 
Lodato and Mechtly (1971). The operating frequency used in the calculations is 
2.66 MHz. Figure 3.2 shows the phase change due to reflection as a function of 
altitude. The variation in the phase difference due to reflection for these col- 
lision frequencies is small. In Figure 3.3 the expected values of A<}> are plotted 
for the same three rocket profiles. Because electron densities are low at lower 
altitudes, the values of A<j> expected at lower altitudes are near zero. At higher 
altitudes the phase due to the difference in refractive indices becomes large. 

In making differential phase measurements, then, one would expect to obtain 

(A4> 0 - A<j>^) versus. altitude profiles similar to those in Figure 3.3. 

3.2 Volume Scattering Theory 

The theory of volume scattering has been proposed as an improved model for 
partial reflections. Flood (1968) first applied volume scattering theory to the 
D-region partial-reflection experiment. He was able to calculate electron con- 
centrations from differential absorption using this theory. For their differ- 
ential phase measurements, von Biel, et_ al_. (1970) used volume scattering the- 
ory for calculating electron concentrations. In addition, the results of their 
measurements were offered as evidence that the volume scattering theory was more 
accurate than the simpler Fresnel theory. 

In contrast to Fresnel theory, the electromagnetic waves incident on the 
ionosphere are backscattered by a volume of random fluctuations in refractive 
index. For experimental calculations the size of the volume is taken to be the 
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Figure 3.1 Electron concentration profiles from rocket measurements (after Mechtly and 
Smith, 1968) 
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PHASE OF REFLECTION COEFFICIENT (DEGREES) 


Figure 3.2 Difference between phase components of the reflection 
coefficients calculated from rocket measurements 
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product of the area defined by the antenna beamwidth and the resolution length 
of the transmitter pulse, cx/2. The important contribution of volume scattering 
is that the changes in refractive index throughout the scattering volume are not 
negligible as assumed in Fresnel theory but affect the absorption and phase of 
the scattered wave. 

In the experiment of von Biel, et_ al_. (1970), an expression for the corre- 
lation of the complex electric field intensities of each magneto-ionic mode was 
given as the starting point. For the purpose of comparing the two theories the 
conjugate of the expression for the correlation coefficient given by von Biel 
will be used here. 


P(h) 


<e E*y 

x o x 

KIe I 2 ><Ie x i 2 » 


1/2 


(3.14) 


where the angle brackets denote ensemble averages. The correlation coefficient 
can be separated into its real and imaginary parts. The ratio of the imaginary 
part to the real part can be taken from von Biel, e£ a]_. (1970) 


pj(h) = 

p T (h) 



(n* - n£)dz + p (h) 


(3.15) 


The term 3(h) includes the effects on radio waves due to backscatter from the 
ionosphere. 3(h) is given by 


3(h) = (y x 


- V 


sin(otA) exp(-ag) 

1 - cos(aA) exp(-ag) 


(3.16) 


where y Q and y^ are the phase of the reflection coefficients defined in (3.3) 


through (3.6) 
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„ i i 
5 = n + n 

o 

x 

a = 2 itct/A 

t = transmitted pulse length in seconds. 

Only the last term of (3.16) is a function of the pulse length x. The limit of 
the last term as x approaches zero can be found by using L' Hospital's rule. 


£im sin(otA) exp(-ag) _ A_ 

x->0 1.0 - cos (aA) exp(-aO _ C 


(3.17) 


Therefore, in the limit as x approaches zero Equation (3.15) becomes 


p 1 (h) = 2u 

r,, . c 

P (h) 


(n r - n r )dz + (y - y ) • 
o x' 'x ' o J 


(3.18) 


The right hand side of this expression is identical to that of Equation (3.9). 

The above derivation is a check on the expressions derived from Fresnel theory 
since the two theories should give the same result as the pulse length becomes 
close to zero. 

In order to implement a differential phase system to calculate electron con- 
centrations a reflection model must be chosen. The electron-density profiles of 
von Biel, et_ al_. (1970) and von Biel (1971) were calculated using volume scatter- 
ing theory. The complexity of the expressions necessitates the use of iterative 
techniques to determine electron concentrations. The Fresnel model will be used 
here because of the simpler expression for electron density. While there may be 
some question as to which theory gives the best results, the Fresnel theory will 
give results sufficiently accurate to evaluate the system. 
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4. DIFFERENTIAL PHASE MEASURING SYSTEM 

A method for measuring the phase difference between the ordinary and extra- 
ordinary modes was developed by Parkinson (1955) and later improved by Ferraro 
(1959) to study the polarization of ionospherically-reflected LF waves. A simi- 
lar method was applied to differential phase partial -reflection measurements by 
von Biel, et^ al^. (1970). Two orthogonal antennas in the x and y directions 
(x x y = z) are fed such that the amplitude and phase of the wave to each antenna 
can be controlled. If the amplitudes are equal and the phase of the wave in one 
antenna leads the other by 90°, the antennas can transmit and receive circularly 
polarized waves. Both right and left circular polarization can be obtained in 
this manner. If the waves are in phase at the antennas with equal amplitudes or 
if only one antenna is fed, the resulting polarization will be linear. The mea- 
sured values are amplitudes. When these amplitudes are combined in the proper 
way, the phase change due to ionization in the D region can be evaluated. 

Parkinson (1955) and Ferraro (1959) used a method of calculation which 
utilizes the angle of phase difference between the ordinary and extraordinary 
waves. On the other hand, von Biel, et al_. (1970) calculated a correlation co- 
efficient. Both methods give the same result, as will be shown; however, each 
method can give different insight into the measurement. 

4. 1 Angle of Phase Difference 

The complex intensity of the backscattered field from a height h is 
t - K * Ci * jy) * (X - jy) 


(4.1) 
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The transmitted wave is linearly polarized and pulsed. Mode coupling and oblique 
echo effects are assumed to be negligible. From (4.1) the electric vector of the 
ordinary ray as a function of time [exp(-ja)t)] is 

A 

E = — (x coswt + y siniot) • (4.2) 

° /J - 


Similarly, 


£ = — [x cos(wt + <j>) - y sin(ot + <£)] (4.3) 

x ^ 

is the electric vector of the extraordinary ray. Both the ordinary and extraor- 
dinary rays experience a phase shift due to ionization and reflection. However, 

since only the total phase difference 4> is needed, (f> is added to the argument 
of the extraordinary mode while the phase of the ordinary wave is referenced to 
zero. The angle cf> is the quantity to be measured. 

The angle <J> is measured by combining the amplitudes of received partial re- 
flections for a suitable set of polarizations. Let 
A q - Ordinary wave transmission and reception 
A^ - Extraordinary wave transmission and reception 
The next four amplitudes are received following a linearly polarized transmission. 
A^ - Linearly polarized reception in x direction 

A^ - Linearly polarized reception in y direction 

A^ - Linearly polarized' reception midway between x and y directions 

A A 

A^ - Linearly polarized reception midway between -x and y directions 
The voltage at the terminals of the receiving antenna oriented in the x 


direction is 
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O X 

V, = cosiot + cos(wt+ <f>) 

1 JT~- ST 


C4.4) 


The amplitude is 


A A 

2 

A 

— + — COS(J> 

+ 

— sind> 

/2 

* 


s? 

\ * 


2 2 
o 
2 


ST— + A A COSd) 

2 o x T 


(4.5) 


Similarly for the receiving antenna oriented in the y direction the amplitude is 

.2 


? A 2 

2 o 


A A cosd) 
ox Y 


(4.6) 


at this point, by measuring A , A , A. , and A~, an expression for the angle <f> can 

O X X ^ 

2 2 

be obtained by taking the difference between A^ and A 2 


cos<|> = (A 


2 .2 


1 2 


a;)/(2a o • A x ) 


(4.7) 


2 2 

The sum of A^ and A ^ is 


,2 A 2 A 2 A 2 

A. + A_ = A + A 
1 2 o x 


(4.8) 


Expressing coscj> as 


cos<j> = 1/2 


r A 2 - A^ 
1 A 2 


A 2 ^ A 2 

L*! * A 2 . 


.-,2 , 2 ' 

A + A 
o x 


, A -A 
’-ox 


(4.9) 
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gives a relation which is independent of system gain changes between circular 
and linear polarizations. 

If the polarization of the receiving antenna is rotated such that it is 
polarized in a direction midway between x and y, the angle <J> is effectively in- 
creased by tt/2. Therefore, the electric vectors are 

I: = A [x coswt + y sin;Dt]//"T 

o o J 

i = A^[x cos (art + 4> + tt/2) - y sin(u)t + </> + 7r/2)]// _ 2" . (4.10) 


The voltage from the x component is given by 


O X 

V, = cosojt + — cos(wt + 4> + tt/2) 

3 yfl /2 


The amplitude is given by 



A A sinb . 
ox 


(4.11) 


(4.12) 


By orienting the polarization midway between -x and y, the angle <j> is effectively 

3tt 

increased by — . The amplitude of the fourth polarization is given by 


A 


2 

4 



A A sin* 
ox 


(4.13) 


Again an expression for the angle <j> can be found by taking the difference between 


A^ and A^ 
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sin<j> 


. 2 .2 

( A 4 ~ A 3) 
2 A • A 


(4.14) 


2 2 

By inserting the sum of A^ and A^ into (4.14) as was done in (4.9) 


sin<j> = 1/2 


r A^ - A 2 
4 3 


.2 A 2 
l A 4 + A 3 


'A 2 A 2 ' 

A + A 
o x 


(4.15) 


Equations (4.9) and (4.15) are sufficient to find <p , provided the amplitudes are 
sampled sufficiently close together. However by taking 


tan* • 2iSt = 
Y COScp 


A 2 A 2 

A 4 " A 3 


A 2 A 2 

A 4 + A 3- 


r A 2 
A 1 


, 2-i 


A 2 . A 2 

A 1 + A 2 


(4.16) 


and retaining the signs of sin<|> and cos<|>, the difficulty of locating the correct 
quadrant is eliminated. 

Von Biel, et al . (1970) have derived (4.9) and (4.15) in terms of the corre- 
lation of the ordinary and extraordinary backscattered fields. The correlation 
coefficient (Equation 3.14) is given by 


P (h) 


<E E 
\ o x / 


«|e 0 i 2 > <iE x r» 


1/2 


(4.17) 


where the angle brackets denote ensemble averages. From Equation (4.1), A^ and 
A 2 can be written as the voltages in the x and y directions, respectively. 
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(4.18) 


Solving for A Q and in terms of A^ and 


A 

o 



j a 2 ) 


A 

x 



j a 2 ) 


(4.19) 


Substituting equations (4.19) into (4.17), 


<v; > 


2~r~i72 


K A „ >•<*,» 


<l4 (*! - )V • 4 (A 1 - i a 2 )] > 

C<A 0 2 > • <A x 2 »V2 


1/2C<A 1 2 > - <A 2 2 > ) - j<Aj A 2 > 

k \ 2 > • <\ 2 » 1/2 


(4.20) 


Amplitudes A^ and A^ can be written from their definitions. 


. _ S2 , /2 . 

A 3 2 A 1 + 2 A 2 


. /2 . S2 . 

A 4 2 A 1 + 2 A 2 


(4.21) 
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Solving for and A ^ and taking the product. 


> A = 
A 2 


1 

n 


(A 3 ♦ 


V * /I (A 3 ' V 




Separating (4.20) into its real and imaginary components, 


r 

P 



i , 1 

P (h) = j 


<\ 2 > - < A 2 2 > 


‘< A o 2 >*<\ 2 > 

,_<\ 2 >* < A 2 2 > j 


_(< A 0 2 >-<\ 2 » 1/2 

"< A 4 2 > - < A 3 2 > 


r <C> * <\ 2 > ’ 

_< A 4 2 > + < A 3 2 > 


_(< A „ 2 >‘<\ 2 » l/2 _ 


again since. 




+ 


2 


(4.22) 


(4.23) 


(4.24) 


the expressions (4.23) have been written as ratios in order to eliminate errors 
due to system gain changes. Equations (4.23) are identical to Equations (4.9) 
and (4.15). 

4. 2 Differential Phase Measuring System 

The partial-reflection differential phase measuring system is diagrammed in 
block form in Figure 4.1. The system is a modified version of the one used by 
Pimat and Bowhill (1968), Reynolds and Sechrist (1970), and Birley and Sechrist 
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Figure 4.1 Block diagram of differential phase partial reflection system 
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(1971). Much of the detail of the equipment such as the receiver and trans- 
mitter has been described by Henry (1966). 

Pulser 

The pulser controls the system by synchronizing the various components of 
the system. Two pulses 33 milliseconds apart are supplied to trigger the trans- 
mitter, control the polarizations of the polarization networks, and initiate 
the computer data collection. The repetition rate for the set of two pulses is 
variable, but is generally at 1.25 sets per second. 

Transmitter 

The transmitter operates at a frequency 2.66 MHz in a pulsed mode. The 
pulse width is continuously variable, but is usually set at 50ys or 25ys. The 
output power of the transmitter is 40 kilowatts. 

Transmitter Antenna Polarization Network 

Figure 4.2 is a block diagram of the network designed to polarize the anten- 
nas. The amplitudes of the polarizations desired--A Q , A , A^, A^, and A^-- 

were described previously and their corresponding polarizations are illustrated 
in Figure 4.. 4. The transmitter power is delivered to a divider (Figure 4.2) 
which sends one half of the signal to each feed line. An attenuator in each line 
equalizes the power that will be delivered to the respective antennas. To obtain 
circular polarization as required for the ordinary (A q ) and extraordinary (A ) 
modes, a pair of crossed dipoles must be fed with waves 90° out of phase with each 
other. A 180° phase shift in an antenna feed line will give circular polarization 
in the opposite direction. If the 90° phase delay is deleted, linear polarization 
results. The direction of linear polarization is controlled manually. With the 
phase delay deleted and power fed to both antennas, the polarization is directed 
northwest as required to obtain A^ and A With the west -east antennas deleted 
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the polarization is directed towards the north, for A^ and A^. 

Antennas 

Two large antenna arrays, one for transmitting and one for receiving, are 
used in the partial -ref lection system. Each antenna consists of a grid of 60 
half-wave dipoles, 30 in the north-south direction and 30 in the east -west 
direction. The gain is directed vertically with a calculated 3dB beamwidth of 
14°. 

Receiver Antenna Polarization Network 

The network used to polarize the receiving antennas is diagrammed in Figure 
4.3. The operation of the receiving network is similar to that of the trans- 
mitter polarization network. A manual phase delete switch has been added to 
obtain linear polarization. Figure 4.4 shows the directions of polarization 
during reception. For the first linear polarization (A^), the direction of 
polarization is north-west. A 180° phase shift of the north-south antenna 
changes the direction to south-west, which gives the required direction for re- 
ception of A. 2* The direction of polarization for A^ and A^ is identical to that 
of A^ and A^, respectively. 

Receiver 

A solid-state, superheterodyne receiver is used in the partial -reflection 
system. The center frequency is at 2.66 MHz with a bandwidth of about 40 MHz. 

The gain is about 80 dB. Provision is made for adjusting the RF and IF gain 
separately. A manually- switched attenuator is used in front of the receiver to 
attenuate signals that would otherwise saturate the receiver. A high-pass filter 
with a cutoff frequency of 1.9 MHz is used to eliminate interference from local 
broadcast stations. An oscilloscope is used to monitor the partial reflections 
at the receiver output. 
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Analog-to-Digital Converter and Computer 

An HP5610A analog-to-digital converter (ADC) digitizes the dc output of 
the receiver. Sampling at a rate of 100 kHz, the ADC has a height resolution 
of 1.5 kilometers. The ADC has a voltage resolution of 9 bits for an input 
voltage range of from 0 to 1.0 volts. An encode pulse generator, which is 
triggered by the pulser, gates the ADC such that it digitizes the voltages for 
the heights of interest. 

A PDP 15/30 small computer is used to store and process data. The computer 
has 16K words of core where each word consists of 18 bits. Data are stored on 
DECtapes and can be processed at a later time. Most of the programming is done 
in FORTRAN with the exception of some special programming which must be done in 
MACRO or assembly language. A more complete description of the computer and its 
peripherals was presented by Birley and Sechrist (1971). 

Data Collection and Processing Programs 

Data are collected from the partial -reflection system through the use of 
the computer program, DLOGD. Each pulse from the pulser generates a frame of 
samples which are digitized by the ADC and stored in memory. Each frame con- 
sists of four noise samples where no partial reflections are expected and 21 
signal samples, covering the altitude range from 60 to 90 kilometers in 1.5 kilo- 
meter increments. Two frames, for example, one for the ordinary mode and one for 
the extraordinary mode, constitute a set. When nine sets are present in memory, 
the data are written out on DECtape for storage and later processing. This 
sequence is repeated until a sufficient number of samples have been collected. 

-Processing is controlled by MASTER. MASTER handles input -output instructions 
from the operator that allow him to select the files to be used and dictate how 
the data are to be processed. MASTER calls the subroutine PR0AXD which does the 
bulk of the processing. PR0AXD reads the values of the data off the DECtape and 



calculates the mean square of the amplitudes. Noise is extracted from the sig- 
nal using the noise samples taken during data collection according to the fol- 
lowing scheme 




S+N 


' 4 


where 

Ag is the amplitude of the signal; 

Ag +N is the amplitude of the signal plus noise; and 
A^ is the amplitude of the noise. 

The average values of the amplitudes for A q , A , A^, A^, A^, and A^ for each 
height are stored in memory. PROAXD, then, calls subroutines to calculate elec- 
tron concentrations from the ratios of A x to A q (see Birley and Sechrist, (1971), 
for a description and listing of the subroutines) . DIFASE is called by MASTER 
to use Equations (4.16) and (4.23) for obtaining the desired phase difference 
angles and correlation coefficients. The electron concentration, as expressed 
in Equation (3.13) is calculated by a straightforward FORTRAN program called 
PHELD. 

Operation of the System 

Data are collected over a 15-minute period. The collection period is di- 
vided into three 5-minute segments. During the first segment the polarization 
shift networks are set for circular polarizations and the amplitudes of the ordi- 
nary and extraordinary modes are collected and stored in one data file. During 
the second 5-minute segment, the 90° phase delays in both polarization networks 
are deleted manually. Two linear polarizations (A^ and A^) are collected and 
written on a second data file. During the final segment, one part of the 



transmitting antenna is switched to a dummy load, changing the direction and 
amplitude of the polarization. The remaining two polarizations (A^ and 
are then recorded on a third data file. Three hundred seventy-eight samples 
of each of the six amplitudes are compiled over the entire collecting period 
The data are then ready for processing. 
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5. RESULTS OF DIFFERENTIAL PHASE MEASUREMENTS 

The objective of making differential phase partial -reflection measurements 
is to obtain electron -concentration profiles of the D region. In this chapter 
several profiles will be presented to demonstrate the results of the technique. 
The data collection period began on October 13, 1971, and continued through 
November 5, 1971. The location of the system is northeast of Urbana, Illinois 
at geographic coordinates: 40°N, 88°W. 

The advantage of making differential phase measurements by the system des- 
cribed in Chapter 4 is that differential absorption measurements are made simul- 
taneously. Comparisons can be made quite readily of the results from the two 
techniques. A common collision-frequency profile was chosen for calculating 
electron concentrations by both methods: this collision-frequency profile (Fig- 

ure 5.1) was chosen as representative of several equinox profiles presented by 
Lodato and Mechtly (1971) . 

The phase difference, the ratio A x /A q , and the amplitude of the correlation 
coefficient are plotted as a function of altitude for two data runs in Figures 
5.2 and 5.3. These phase difference profiles have the same general shape as 
those presented in Figure 3.3. The phase difference increases slowly at the 
lower altitudes, whereas at higher altitudes it increases very rapidly. The 
rapid increases in phase difference occur at lower altitudes for the data mea- 
sured by the system than for the calculated values from the rocket data. 

The amplitude of the correlation coefficient in Figures 5.2 and 5.3 is a 
decreasing function of altitude. This indicates that the ordinary and extra- 
ordinary waves tend to become uncorrelated as altitude increases. Ionospheric 
characteristics such as electron concentrations and the nature of the partial 
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reflections can be deduced from the correlation coefficient. From the results of 
von Biel, £t^ al_. (1970), one would expect the correlation coefficient to be near 
unity at low altitudes. The measured values of p were seldom larger than 0.8 
below 75 kilometers. Since p is a function of the refractive indices of the waves, 
it is a function of the electron concentration. If the values of p are too small, 
the electron concentration from p would be too large. The correlation coefficient 
shows a rapid change at high altitudes as do the phase difference and A^/A q . The 
use of the amplitude of the correlation coefficient to calculate electron concen- 
trations was suggested by von Biel (1971). 

The electron-density profiles from differential phase and differential ab- 
sorption are presented in Figures 5.4 through' 5.10. In general, the agreement 
between differential phase and differential absorption results is quite good. 

Values of electron concentration were obtained over altitudes from 70 kilometers 
to 86 kilometers depending on the kind of day. In particular the profiles from 
absorption and phase of Figure 5.7 show good agreement over a range of 10 kilo- 
meters. The solar zenith angles were less than 60° for all of the measurements. 

The length of the transmitted pusle was either 50ys or 25ys. The first runs 
of data were collected with a 50ys pulse. Later 25ys was tried and seemed to give 
somewhat better results. 

The seven profiles show some variability in electron concentration over the 
data collection period. The values of the electron density are large at low 
altitudes during mid-October. Beginning with Figure 5.8 the electron density 
appears to have decreased. These observations are consistent with those made 
by the regular differential absorption measurements during this period. 
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6. DIFFERENTIAL PHASE SYSTEM EVALUATION 

A differential phase partial -reflection system has been developed, and the 
results from it have been presented in Chapter 5. It is necessary to examine 
the system and the results to arrive at some quantitative conclusions as to the 
merits of the system. To evaluate the system, sources of error must be consid- 
ered. The magnitude of these errors and their cumulative effect on the final 
results must be analyzed. Suggestions can then be made for minimizing the 
errors and improving the overall quality of the system. From this evaluation 
one can better understand the merits and limitations of the system. 

The sources of error can be divided into two categories : those errors due 

to uncertainties in measuring techniques, and those due to faulty assumptions in 
the theoretical calculations . Several sources of error which apply to this 
system will be considered under each of these categories. 

One of the best ways of evaluating a differential phase system is to compare 
the results with those from differential absorption measurements made simulta- 
neously. This comparison has been done in Chapter 5. Further comparisons with 
data from other types of experiments would also be worthwhile and will be done 
here. 

6.1 Uncertainties in Measurement Techniques 

In Chapter 4 the polarizations required for each of the six amplitudes were 
assumed to be exactly circular or exactly linear as required. Whether or not 
the polarizations are really circular or linear must be considered as well as 
what effect the non- circular or non-linear polarizations will have on the mea- 
surement. The criteria for transmission or reception of circular polarization 
are that the amount of attenuation of the wave in each feedline be equal and 
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the wave of one feedline be 90° out of phase with respect to the wave in the other 
feedline. For linear polarization, the attenuation of each wave must be the 
same, and there must be zero phase difference between the two waves. If the 
wave attenuation in each feedline is different, the angle <j> in the equations for 
each amplitude (4.5, 4.6, 4.12, and 4.13) has an error. Since the wave attenua- 
tion in each line can be adjusted to within 3%, the maximum possible error in 4> 
is 2°. The maximum error occurs when <j> is an odd multiple of tt/2. As <j> approach- 
es zero the error in <f> becomes nearly zero. Therefore, the effect of an error 
in the direction of polarization will not be important. A method of testing the 
phase between the two feedlines is to measure the rejection of one mode while 
transmitting the other mode. Each antenna had a rejection of 20 dB for both 
modes. If for linear polarizations the phase between the two feedlines is not 
zero, the polarizations will be elliptical. Again the possible error is less 
than 2°, if the phase difference between the feedlines is less than one degree. 

The output of the receiver is non-linear. If this non-linearity remained 
in the data, the changes in <j> over height would be significantly affected. A 
computer program which approximates the receiver output by a piece-wise linear 
curve compensates for the non-linearity. 

When the gain of the receiver is set for measuring partial reflections at 
low altitudes, the signals from the stronger reflections at higher altitudes 
saturate the receiver. If the data from saturated signals are accepted, the 
mean squared amplitude will be too low. Therefore, samples that saturate. the . 
receiver must be rejected. PROAXD rejects saturated signals and lists the 
number of rejections for each amplitude. If the number of rejections at a 
particular altitude is significant, that is, on the order of 10% of the sam- 
ples, the data are considered to be invalid. 
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The timing of the control pulses in the system is very important for the 
precise determination of altitude. The center of the transmitted pulse is 
taken as the zero height. Signals from 60 kilometers and 90 kilometers are 
received 400 and 600 ys later, respectively. An error of 6.67 ys will result 
in a height error of one kilometer. The present system can be adjusted such 
that the height error is no more than two kilometers. 

The noise which can contaminate the data can be either internal or exter- 
nal. The internal noise in the system is very low. External noise such as 
atmospherics or interference can be very significant. When a received signal 
is contaminated by a large source such as a flash of lightning or a strong 
radio transmitter the entire frame is rejected. By sampling at about 30 kilo- 
meters the amplitude of the noise is recorded. The signal can be subtracted 
from the signal plus noise as described in Chapter 4. The signal -to-noise 
ratio was examined for each height, and data with ratios of less then 1,5 were 
not accepted. 

Each of the six amplitudes measured at a particular height is a random 
number the squared value of which varies about a most probable or mean squared 
value. If the number of samples over time is large, the amplitudes squared 
will tend to a normal distribution about the mean squared value. The standard 
deviation can be estimated by (Carpenter and Bowhill, 1971) 

/c<|a | 2 » 2 (<|A n l 2 >) 2 V /2 

( 5 ♦— 4S ) (6 - 

2 2 
where ^ | | is the mean squared amplitude of the signal, | | y is the 

mean squared amplitude of the noise, and is the number of samples. For the 

data of Figure 5.5, the standard deviation is +_ 5% for each amplitude. 
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The effect of a standard deviation on the phase difference and eventually 

the electron concentration will now be considered. Each of the four mean squared 

2 

amplitudes in Equation (4.16) <^|A^| y can be represented by a number with a 

standard deviation u^ . We will assume that the squares of the four standard de- 
viations are equal . 



The standard deviation of A4> for Equation (4.16) is given by 


( 6 . 2 ) 


where 


o(Acf>) 


^~2 °(P) 

1 + p 



The standard deviation of p is given by 


(6.3) 



(6.4) 
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The standard deviation for the phase difference between two heights is given by 


a(A <|) 2 - A^) 


/ 2 

° (Po) 


° 2 (pJ 


(6.5) 


The standard deviation of the electron density is given by 



( 6 . 6 ) 

These standard deviations for a u of 10% have been calculated from the data of 
October 19, 1971, at x = 55° and are tabulated in Table 6.1. The standard devia- 
tions of the electron concentrations a(N) range from 6% at 75.75 kilometers to 
24% at 80.25 kilometers. 

6.2 Assumptions and Theories 

In order to measure parameters of the ionosphere, one must have some knowledge 
of what the ionosphere is really like. Assumptions are made concerning the ioni- 
zation, the collision frequency and the reflection model. If these assumptions 
are wrong the analysis and results will be wrong. A detailed development of these 
assumptions will not be attempted here, but rather a summary of the problems in- 
volved will be presented. 

The best way to evaluate a system for measuring electron concentrations is 
to compare the results from the system with a real electron density profile. 
However, it is very difficult to know what the real electron density is because 
so few of the parameters necessary to calculate the ionization such as concen- 
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Table 6.1 

Standard deviations calculated from data of October 19, 1971, x = 55°. 


Altitude 

(km) 

a 

(P) 

o(A<j>) 

(Degrees) 

Altitude 

(km) 

ct Ca<j> 2 -A4> ) 
(Degrees) 

a(N). 
(cni ) 

75.0 

1.85 

x 10" 2 

0.466 

75.75 

1.84 

29.2 

76.5 

2.78 

x 10" 1 

1.87 

77.25 

1.80 

25.6 

78.0 

2.56 

x 10 _1 

0.248 

78.75 

5.43 

72.5 

79.5 

1.59 


5.43 

80.25 

12.6 

154. 

81.0 

3.46 

x 10" 1 

11.4 

81.75 

29.6 

359. 


82.5 


9.52 


27.2 
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trations of ionized and neutral species and reaction rates are known at the time 
of measurement. The rocket profiles are considered to be close to the real pro- 
files; however, the number of profiles from rockets are few and the variability 
is great. For example, rocket profiles from 1964 through 1969 (Mechtly, 1972) 

_3 

show a range of 100 to 900 electrons cm at 80 kilometers. All but two of the 
20 profiles from differential phase measurements have values of electron concen- 
tration within this range at 80 kilometers. 

Another factor is the change in ionization over the sampling period. This is 
particularly important in the case of the system used here because amplitudes A^ 
and A ^ are taken 5 minutes apart from amplitudes A^ and A^. The change in ioni- 
zation is largely dependent on the change in solar zenith angle. By minimizing 
the change in solar zenith angle, changes in ionization can be minimized. All 
of the data collected were for solar zenith angles less than or equal to 60°. 
Therefore, changes in solar zenith angle were on the order of one degree or less 
over the 15-minute collection period. If the six amplitudes could be taken simul- 
taneously at a high pulse repetition rate, reliable measurements can be made dur- 
ing periods when the zenith angle is changing rapidly. 

The choice of the proper theory has been considered in Chapter 3. While the 
two methods should give similar results, the difference in the results from the 
two theories may be important. The concept of the correlation coefficient intro- 
duced by von Biel et^ al_. (1970) may prove useful for evaluating the two theories. 
The correlation coefficient is a measure of the relative influence of the medium 
on the two waves. One would expect the two waves to become uncorrelated with 
increasing altitude. The phase of the reflection coefficient is essentially 
the same as the phase difference derived from Fresnel theory. However, the 
magnitude of the correlation coefficient could provide useful information. Fur- 
ther measurements and study of the correlation coefficient are recommended. 
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6.3 Collision Frequency 

The calculation of electron concentrations using Equation (3.12) depends on 
an assumed collision-frequency profile. The profile of Figure 5.1 was used to 
obtain the results presented in Chapter 5 . One explanation for the discrepancy 
between the profiles from differential absorption and those for differential phase 
is that the assumed collision-frequency profile was incorrect. 

The effect of changing the collision frequency can be demonstrated by calcu- 
lating electron concentrations from the differential absorption and the differential 
phase expressions. A^/A^ and A<j> are kept constant as the collision frequency is 
varied. As an example, when the collision frequencies of Figure 5.1 were doubled, 
the electron concentrations from differential absorption decreased by as much as 
75% of their former value. From differential phase the electron concentrations 
increased by as much as 200%. The changes in electron concentration were greater 
for differential absorption for altitudes above about 80 kilometers. For differ- 
ential phase the changes were greater below 80 kilometers. 

Fejer (1961) suggested the use of differential absorption and differential 
phase partial-reflection measurements to calculate collision frequencies. With 
the information presented above we are able to develop a simple algorithm that 
will allow us to deduce a collision-frequency profile. This method assumes that 
given the real collision frequency, the electron concentrations calculated from 
differential phase and differential absorption are equal. Using the profile of 
Figure 5.1 as an initial set of collision frequencies the electron concentrations 
are calculated by both methods. From these electron concentrations a multiplier 
is calculated as follows: 


M = N DA /N Dp 
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where and N^p are the electron concentration from differential absorption and 
differential phase respectively. A new value for the collision frequency is found 
by multiplying the old collision frequency by M. This iterative process is con- 
tinued until the two calculated electron concentrations are within 0.1 electron 

_3 

cm of each other. Approximately six iterations are needed. If this technique 
is employed for several heights, a set of collision frequencies can be plotted. 

A least square fit of this set of points gives the new collision-frequency pro- 
file. 

A collision-frequency profile has been deduced from the data of October 19, 
1971 (Figure 5.7), This profile is the dashed line shown in Figure 6.1. The 
new profile, is a least square fit of the collision frequencies plotted nearby. 

The solid line is the initial profile from Figure 5.1. The new profile is slightly 
higher and is rotated from the original. 
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•:7. CONCLUSIONS AND SUGGESTIONS FOR FURTHER WORK 

A successful demonstration of a differential phase partial -reflection sys- 
tem has been made. Electron -density profiles can be obtained from differential 
phase measurements over a height range of approximately 75 to 85 kilometers. 

The Fresnel theory is sufficient to calculate electron concentrations and to 
evaluate the system. An error analysis is useful for evaluating the system 
and the results obtained from the system. Differential absorption and differ- 
ential phase results can be combined to obtain a collision -frequency profile. 
However, the collision frequency obtained in this manner should be regarded 
with caution. By comparing the deduced collision -frequency profile with ex- 
pected profiles, observations as to the validity of the results can be made. 

The addition of a differential phase measuring system to the existing par- 
tial-reflection experiment would be useful for the study of the lower ionosphere. 
Differential phase results provide a check on the reliability of the differen- 
tial absorption results. A partial -reflection system consisting of both dif- 
ferential absorption and differential phase would be useful for studying both 
D-region ionization and the reflection processes. 

Several improvements of the partial -reflection system in general and dif- 
ferential phase measurements in particular are recommended. An improved receiver 
design would allow more versatility in programming. With a receiver having a 
linear output, the time required for collecting and processing data would be 
decreased by more than 50%. The six pulses should be measured sequentially in- 
stead of in three segments. The present pulser should be replaced by a digital 
input -output register that would, allow the computer to control the system. The 
different amplitudes would be obtained by writing the proper software. 
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Further investigation of the use of volume scattering theory for evaluation 
of differential phase measurements is recommended. The correlation coefficients 
from an assumed electron-density profile should be compared with the phase dif- 
ference calculations. Also, the magnitude of the correlation coefficient should 
be used in conjunction with the phase difference and A^/A^ ratio measurements 
for investigating the performance of the system. 

The use of statistical data analysis for evaluating partial-reflection data 
is recommended. The standard deviations, developed in Chapter 6, give some mea- 
sure of the confidence that can be placed in the measurement. It would be desir- 
able to incorporate this analysis as a routine part of the data processing. 


I 
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Ao 

Ax = 

A1 

A2 

A3 

A4 

HT 

RHD = 
RHOI = 
RHOR = 
PHI = 
DPHI = 
GNU = 
PL 

QL = 
CMEG = 
CTO = 

CTX = 

CFO = 
CFX = 


APPENDIX 

List of symbols used in the computer programs 

amplitude of ordinary polarized waves 
amplitude of extraordinary polarized waves 
linear polarization corresponding to A^ in the text 

linear polarization corresponding to k^ in the text 

linear polarization corresponding to A^ in the text 

linear polarization corresponding to A^ in the text 

altitude in kilometers 

amplitude of the correlation coefficient 
imaginary part of the correlation coefficient 
real part of the correlation coefficient 
A<f> 

A<f>2 " A<f>^ 

v 

m 

U) + OJ^ 

“ - w L 

angular operating frequency 
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TITLE: DIFASE 

********** 
°UR POSE: SUBROUTINE TO CALCULATE DIFFERENTIAL PHASE FROM AO, AX, 

A. 1 , A2 , A3 , AND A4. 

THE ?POPFR QUADRANT IS FOUND BY EXAMINING THE SIGNS OF THE NUMERATOR 
AND DENOMINATOR BEFORE TAKING THE ARCTANGENT. THE CHANGE IN 
DIFFERENTIAL PHASE BETWEEN TWO HEIGHTS IS FOUND BY SUBTRACTING 
THE PHASE AT A LOWER HEIGHT FROM THAT AT A HIGHER HEIGHT. THE 
ANGLES ARE IN DEGREES. 

* **.*.* * * * * * 
SUBROUTI NE DI FA SE (HT ,MM, NN , A 0 ,A X, A 1 ,A2 , A3 ,A 4 ,PHI ) 

DIMENSION HT (5 0) ,AX(50) ,A 0 (50) ,A ! (50),A2(50) ,A3 (50) ,A4 (50) , 

IRHOR (50) ,R HOI (50) ,RHO(50) ,PHI (50) 

DO 25 J=MM,NN 
CIR=A0(J)**2+AX(J>**2 
BLI NR=A1 (J)**2+A2 (J)**2 
BLI Nl =A3 (J)**2+A4 (J) **2 
25 WRITE (6,30) HT (J) ,CIR ,8LI NR # BLI NI 

3 0 FORMA T (5X,F5 .2 ,4X,F 1 0.2 ,4X,F 1 0.2,4X,F 1 0.2) 

WRITE (6,31) 

31 FORMA T ( 1 HI , 1 X, 6HHEI GHT ,4X , RHRHO REAL , 6X, 8HRH0 IMAG,9X,3HRH0, 

1 8X,3HPHI ) 

DO 100 KrMM,MN 
PROD=AO(K)*AX(K) 

R NUM=(A.l (K)**2-A2 (K) **2) *(A0(K)**2+AX(K)**2) 

RDENr (A 1 (K)**2+Al (K)**2)*PR0D 
R HOR (K) =R NUM/(2 .0*RDEN) 

G NUM= (A 4 (K)**2-A3 (K) **2)* (A0(K)**2+A X(K)**2) 

GD E N = ( A 4 ( K ) **2+A 3 (K) **2)*PR OD 
R HOI (K)=GNUM/(2.0*GDFN) 

ARHOI =ABS(RHOI (K) ) 

ARH0RrABS(RH0R (K)) 

IF(RHOI (K) .LT.0.0) GO TO 40 
I F (RHOR (K) ,LT . 0. 0) GO TO 35 
PHI ( K) =5 7.295 78*ATA N(ARHOI /ARHOR) 

GO TO 50 

35 PHI (K) = 1 80.0-5 7 .295 78* A TA N (AR HOI /ARHOR) 

GO TO 50 

40 IF(RHOR(K) .LT.0.0) GO TO 45 

PHI (K) = 3 60, 0-5 7.295 78*ATA N (ARHOI /ARHOR ) 

GO TO 50 

45 PHI (K)= 1 80.0+57,295 78*ATA N(ARHOI /ARHOR) 

5 0 CONTINUE 

RHO(K)=SQRT(RHOI (K)**2+RH0R (K)**2) 

100 WR I TE ( 6 , 75 ) HT ( K ) ,R H OR ( K ) ,RHOI (K) ,RHO(K) ,PHI (K) 

75 FORMAT (2X,F5 .2 ,2X,E 1 2. 8,2X,E 12 . 8,2X,E12 .8 ,2X,F 1 0.4) 

RETUR N 
END 
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TITLE: PROA XD 

* * Jfc * 4c 4c 4c 4c * * 

D !IPPO?E: SUP POUT I ME TO READ DATA FROM DECTAPE, AVERAGE IT, A PD 
CALCULATE ELECTRON CONCENTRATIONS. 

SIP ° PVT I NFS CALLED: D R EAD , CALC 

THE DATA FROM THE DECTAPE ARE SQUARED AND AVERAGED. THE MEAN 
SQUARED NOISE OF EACH POLARIZATION IS SUBTRACTED FROM THE MEAN 
SQUARED AMPLITUDE. THE RMS AVERAGES FOR THE NOISE AMD THE AMPLI- 
TUDES ARE TABULATED FOR EACH ALTITUDE. ELECTRON CONCENTRATIONS 
OOP AX AND AO DATA ONLY ARE CALCULATED BY SUBROUTINE CALC. DREAD 
AND CALC ARE SUBROUTINES USED IN THE AERONOMY LAB FOR DIFFERENTIAL 
ABSORPTION WORK. 

* **4c**4c4c** 

SUBPOUTI NE PROAXD ( FNA M,M ,B MXNS ,R E 1 ,R E2 ,HT ,A VA 0 ,A VA X) 

DIME NSI ON ENAM(P) ,AA (PI) ,AB(21 ) ,AVA0(21 ) ,AVAX(2I ) , 

I XO (2 1 ) , IP JO (2 I ) ,1 R JX (2 1 ) ,DI FNO(4) ,DIFNX(4) ,HT(2I ) , 

200(21) ,BX(2 1 ) ,B NO (4) ,B NX (4) ,B NXNC4) ,AXN(21 ) ,PBN0(4) ,BBNX(4) 

S NO = 0 . 

S NX - 0 • 

I P - 0 
IPNOrO 
IP NX -0 
OMOrO. 

OMX=0. 

DO IS 1 = 1 ,21 
X0CI)=0. 

AVAO(I)=PI. 

AVAX(I)=0. 

A A ( I ) = 0 . 

AB(I)=Q. 

IPJO(I)=G 


IS 

IP JX (I ) = 0 


DO 17 1 = 1 , 


0 NO (I ) = 0 . 

1 2 

B NX ( I ) = 0 • 


CALL VALUE 


CALL DINIT 


CALL ESTAT(2,FNAM,L0G) 

I E (LOG . NE • 0) GO TO 40 
WRITE (6, 35) FNAM 

35 F OR MA T ( SH FILE ,2A5,19H NOT FOUND DM DAT 2) 

RETURN 

40 CALL SEEK (2,FNA M) 

KE0E0=0 . 

KEOFX =0 
ID=Q 

4 S DO 1 4 P 1 = 1 ,21 

B 0 (I ) = AA (I ) 

1 4 P BX( I)=A3(I) 

DO 149 1=1,4 
BP NO (I ) =B NO (I ) 
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149 BP.NX(I)=BNX(I) 

CALL DREAD (A A ,9 NO , I ERR ,ID,KEOFO) >,j 

IFCKEOFO.EO.l ) GO TO 50 / 

CALL DREAD CAB ,P NX , I ERR »I D ,KEOFX) / 

I F (KFOFX .FO , 1 ) GO TO 49 / 

601 RMEANO=0. J 

BMEANX=0. • 

DO 440 1=1,4 
P MEA NO =B MEA MO+B NO ( I ) **2 
440 P MEA MX =3 MEA NX+B NX (I ) **2 

BMEA NO=SQRT (BMEA NO/4,) 

BMEA MX =SQRT (3 MEA NX/4 .) 

IFCBMEANO.GT.BMXNS.AMD .BMEA NX.GT .3MXNS) GO TO 510= 

DO 220 1 = 1 ,4 

DI FNO (I ) =B NO (I ) -3BM0 (I ) 

220 DIFNX(I)=BNXCI)-BBNX(I) 

DO 221 1=1 ,2 

J = I+1 
K = I+2 

IFCDIFNO(I) .GT.REl .OR.DI FNOCJ) .GT .RE 1 .OR .DI FNO (K) .GT.RE1 .OR. 
1BM0(J) .GT.510..OR.BNXCJ) .GT.510.) GO TO 71 0 
BMO=BMO+BNO(J)**2 
GO TO 715 

710 I R NO = 1 R N0+ 1 

IF(BNO(J). GT.510.. OR. BNX(J).GT. 510. ) GO TO 720 
715 I FCDIFNXCI) .GT.REl .OR.DI FNXCJ) .GT.RE1 .OR.DI FNX(K) .GT .RE 1 ) 

1 GO TO 72 0 
B MX =8 MX+B NX ( J) **2 
GO TO 221 

720 I R NX =1 R NX+ 1 

221 CONTINUE 

700 DO 47 I =1 ,21 

J=I + 1 
L =1 -1 

I F (L.EQ .0)L= 1 
I F (J.GT .2 1 ) J=2 1 
DIFE=AA (I)-BO(I) 

DI FEA=AA (L)-BO(L) 

D I FEB =A A ( J) -B 0 (J) 

I F (I .GT.9) GO TO 001 

I F (DI FE. GT . RE I. OR. DIFEA. GT.REl. OR.DI FEB. GT.REl .OR. 

1 AA (I ) .GT .51 0. . OR . AB (I ) .GT.510.) GO TO 45 
GO TO 110 

001 IF(DIFE.GT.PE2 .OR .DI FEA . GT ,PE2 . OR .DI FEB . GT .RE2 .OR . 

1AA(I). GT.510.. OR. ABCD.GT. 510. ) GO TO 45 

110 A VA 0 Cl ) =A VA 0 (I ) +AA C I ) ***2 
GO TO 111 

45 IR JO (I ) rlRJO <I) + 1 

IFCAACI) .GT.510. .OR. AB(I) .GT.510. ) GO TO 46 

111 DIFE=AB(I)-BX(I) 

DIFEA=AB(L)-BXCL) 

D I FEB =AB(J)-BX(J) 
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373 

120 

46 

4 7 

C J f? 

c^O (51 


49 
5 °l 


P’3‘4 


5?. 


163 


5 4 


R <7 
<? 7 


IF(I.RT.9.AMD.I.LT.16)00 TO 303 

IF(DIFE.GT.RE1 .OP.DIFEA.fiT.REl .OP.DIFED.GT.REDGO TO 46 
GO TO 120 

IF(DIFE.GT.PE2.0R.DIFEA .GT .RE2 .OR .DIFEB.GT.RE2) 00 TO 46 
A VA X ( I ) =A VA X ( I ) +AB ( I ) **2 
no TO 47 

IRJX(I)=IPJX(I)+1 
CONTI NUE 
no TO 520 
IP rIR+I 

S NO =D MEA M0**2+S MO 
S NX =BMEA NX**2+S NX 
RO TO 4F 
I D =ID - 1 
10 - 10 / 2 . 

P.ID=ID 

DO *3 0 I r 1 , 2 1 

IRJO(I) = IPJO(I)+IR 

IRJX(I)=IRJX(I)+IP 

IR.NC=IPN0+2*IR 

IPNX=IRNX+2*IP 

A V/ MO r SORT ( S NO/BID ) 

AVNX=SQRT(SNX/BID) 

R NO =1 D *2-1 R MO 
P MX r I D*2-I P MX 
DO 52 1=1 ,21 
RSAMO = ID-IRJO(I> 

RSAWX=ID-IRJX(I> 

A VOC =A VA 0 ( I ) /R SAXO-8 MO /P MO 
A VXC =A VA X (I ) /R SA MX-3 MX/R MX 
AVAO(I) = (ABSCAVOC) /A VOC)*SORT (ABSCAVOC ) ) 

A VAX(I ) = (ABS (A VXC) /AVXC)*SQRT (ABS(AVXC) ) 

CONTINUE 

PX0=SORT (BMO/P NO) 

PMX^SORT (BMX/R NX) 

CALL HEAD ( 1 ) 

VR I TE ( 6 , 1 53 ) A VMO ,A VMX,RMO ,1 R NO ,BMX,I R MX 

FOPMA T (2X, 1 7HA VERARE NOISE ORD ,FF. 1 /2X, 1 7HA VERA RE NOISE EXT , 1 X 

1 ,F7.I //2X,9H NOISE ORD ,F F. 1 ,1 5 , 1 1 H REJECTI 0NS/2X, 

29KN0ISE EXT, F«.l ,15,1 1H REJECTIONS) 

WR I TF C 6, 5 4 ) I D , I P 

FOR MA T ( //I X, I 4 , l 4H SAMPLES TAKEN ,5X, I 5 , 1 6H FPAMES REJECTED// 
19K REJECTED ,2X,9K REJECTED ,2X,6HHEI RHT,2X,6HA V, AO, 2 X, 

2 6HA V . A X/4X , 3 H OP D , PX , 3 HEXT ) 

HTI N = 5 P . 5 

DO 53 1 = 1 ,21 
XZ = I 

HT ( I ) = HTI N+XZ* 1 .5 

WPITE(6,5P)IRJ0(I) ,IRJX(I) ,HT(I) ,AVAO(I) ,AVAXCI) 

FORMAT ( 3 X , I 4,7X,I 4,3X,F5« 1 , 3 X , F 6 . 1 ,2X,F6.1 ) 

CONTINUE 

I F ( M.EO .2) RO TO 62 
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DO 5 9 1 = 1,21 

IF(AVACKI) .LE.O.O.OR.AVAX(I) .LE. 0.8)00 TO 59 
XO(I)=AVAX(I) /AVAO(I) 

S9 CONTINUE 

CALL CALC(XO,l ,20,RESP) 

G2. R etur r 

END 


TITLE: f^STER 

********* 

°U R POSE: DATA PP OPES SI NO CONTROL PROGRAM 
SUBROUTINES CALLED: PR OA XD ,D I FA SE 

PROCESSING INFORMATION SUCH AS THE DATAFILE AND NOISE CRITERIA 
IS PEAD INTO THE PROGRAM. PROAXD IS CALLED 3 TIMES. THE FIRST 
TIME AX AND AO DATA APE PROCESSED AND AN ELECTRON CONCENTRATION 
PROFILE IS CALCULATED. THE SECOND TIME A1 AND A2 DATA ARE PRO- 
CESSED, TABULATED, AND STORED IN AN ARRAY. THE THIRD TIME 
A3 AND A 4 DATA ARE PROCESSED, TABULATED, AND STORED IN AN ARRAY. 
DIEASE IS THEN CALLED TO CALCULATE DIFFERENTIAL PHASE FROM AO, 

AX ,A 1 , A2 ,A 3 , AND A4. 

********* 
DIMENSION FNAMC2) ,A A (21 ) ,AB (21 ) ,A VA 0(21 ) ,AVAX(21) , 

1X0(21 ) ,IRJ0(2! ) ,IRJX(21) ,DIFN0(4) ,PIFMX(4) ,B 0 (21 ) ,9 X(2 1 ) , 

2B NO (4 ) ,B NX (4 ) ,P NXN (4 ) , A XN(21 ) ,B BNO (4 ) ,B B NX (4 ) ,A0(21) , 

3 A X (2 1 > ,A 1 (21 ) ,A2 (21 ) ,A3 (21) ,A4 (21 ) ,HT (21 ) ,R HOR (2 1 ) ,R HOI ( 2 1 ) , 
4RH0 (2 1 ) ,PHI (21) ,D PHI (21) 

5 CALL HEAD (0) 

WRITE (5,10) 

10 FORMAT (15H WHICH DATAFILE) 

READ (4 ,2 0) FNAM 
20 FORMAT (2A5) 

Mzl 

WRI TF( 5,25) 

25 FOR MA T (24H COLLISION PR EC . PROPILE/ 

1 25H SUMMER, WINTER, OP EPUINOX) 

READ (4,30)PESP 
30 FORMAT (A5) 

WRITE (5,35) 

35 FORMA T ( 1 4H MA XI MUM NOISE) 

READ ( 4 , 4 0 ) B MX NS 
40 POP MAT (FI 0.0) 

WRITE (5,45) 

45 FORMAT (2 1H REJECTION BELOW 72XM) 

READ (4, 40) RE 1 
WRITE (6,5 0) 

5 0 FORMAT (21 H REJECTION ABOVE 72KM) 

PEAD (4 ,4 0)RP2 

CALL PROAXD (FNAM, M,PMXNS,RE1 ,F E2 , H T , A 0 , A, X ) 

M = 2 

WR I TE ( 6 , 1 0) 

PEAD (4,20) FNAM 
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CALL PROAXD (FNAM,M,RMXNS,RE1 ,RF2 ,HT,A 1 ,A2) 
WRITE (6, 1 0) 

READ (A ,2 PI) FMAM 

CALL PROAXD (FNAM,M,BMXNS,RE1 ,RE2 ,HT,A3,A4) 
WR I TECS, 55) 

55 F OR MA T ( 1 5 H INITIAL HEIGHT) 

READ(4,40)HT! 

WRITF(6,60) 

G0 FORMA T(13H FINAL HEIGHT) 

R EAD (4 , 4 0) HT2 
MM=(KTl-5*.5)/l .5 
NfJr (HT2-5F.5) /I .5 

CALL DIFASE(HT,MM,NN,A0,AX,A1 ,A2,A3 , A 4, PHI ) 
WRITE (6,65) 

G5 FORMAT ( ///7H HEIGHT,2X,16HPHASE DIFFERENCE) 

l=nn-i 

DO 70 LM=MM,L 
MK=LM+1 

DPHI (LM)=PHI (MK)-PHI (LM) 

HTD=HT(LM)+0.75 

70 WFITE(6, 75)HTD,DPHI (LM) 

7G F0RMAT(2X,F5.2,4X,F 10.4) 

CO TO 5- 
F ND 


C TITLE: PHELD 

C * * * * * * * * * 

C PURPOSE: CALCULATE ELECTRON CONCENTRATIONS FROM PHASE DATA 
CAN ASSUMED COLLISION FREQUENCY PROFILE IS STORED IN MEMORY. THE 
C PHASE DIFFERENCE IS READ IN SEQUENTIALLY FROM THE LOWEST HEIGHT. 

C THE QUASI -LONGITUDI NAL APPROXIMATION TO THE GENERALIZED FORMULA 
C FOR THE REFRACTIVE INDEX IS USED. THE FRESNEL MODEL WAS USED TO 
C CALCULATE REFLECTION COEFFICIENTS. THE ELECTRON CONCENTRATION IS 
C PRINTED OUT IN ELECTRONS PER CUBIC CENTIMETER. 

C ********* 

DIMENSION DPHI (20) ,GNU(21 > ,EN(20> 

C3 (X)=(X*(X*(X*(X+24.6531 15)+! 13.54165+1 1 .2*75 1 3 ) + . 0239*34 74 ) 

1 /(X* (X* (X* (X* (X* (X+24 . 65 6*19 )+ 120. 49512) +2 *9. 5*0 *5 )+I49. 2 1254) 
2+9. 3*7273 72)+. 01 *0641 2R> 

C5(X)=(X*(X*(X+6. 6945939)+! 6.901002)+! .1630641 ) /(X* (X* (X* (X* (X+ 
16. 63144 97 ) +35. 355 25 7)+6*. 92 05 05 )+64 . 093464)+4 .3605 732) 

C ENTER COLLISION FREQUENCY PROFILE 
GNU(1)=2.20E7 
G NU (2 ) = 1 . 72 E7 
GNU(3) = 1 .32E7 
GNU(4) = 1 .04E7 
GNU (5) =8.0E6 
GNU (6) =6.2 E6 
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G NU ( 7) =4 . 8E6 
GNU(R)=3.7EG 
GNU (9)=2 .85E6 
GNU (10) =2 .2 EG 
GNU ( 1 1 ) = 1 .70E6 
G NU ( 1 2 ) - 1 • 3 EG 
GNU (13) = 1 .0EG 
GNU(J4) = 7, PE5 
G NU ( 1 5) =G .0E5 
GNU(16)=4.SE5 
G NU ( 1 7) =3 .6E5 
GNU (18) =2. 75 E5 
GNU ( 1 9) =2 . 1 2 E5 
GNU (20) = 1 .G5E5 
G NU (2 1 ) = 1 .35E5 
C ENTER CONSTANTS 

PL=2.59G1 4E+ 7 
QL=7.3886E+6 
E=l .G0210E-19 
EM= 9. 1 090PE-3 1 
PERM=g.854E-12 
0MEG=2. 0*3. 14159*2. 66E+6 
10 WRITE (G, 100) 

100 FORMAT ( 1 HI , 1 5H INITIAL HEIGHT) 

P EAD (4,1 1 0)HT 1 
110 FORMAT(FI0.0) 

WRITE (6,120) 

12 0 FORMAT (13H FINAL HEIGHT) 

READ (4,1 1 0)HT2 
M=(HTl-58.5)/1.5 
N= (HT2-5 8.5 ) /I .5 
T =M- 1 

HT =59 ,25+T* 1 .5 
DO 20 L=M,N 

2 0 READ (4,25) DPHI (L ) 

25 FORMAT(F10.G) 

WR I TECG, 35 ) 

35 FORMAT(lHl) 

WR I TE ( G , 40) 

40 F0RMAT(7H HEIGHT, 4X,9HPHASE DI F,5X, 9HC0LL FREQ , 4X , 9HELEC CONC) 

DO 50 J=M,N 
HT =HT+ 1 .5 

C SEN-WYLLER EQUATIONS 

GNU1 = (GNU(J) + GNU(J+1 ))/2.0 
0 1 =PL/GNU (J) 

02 = PL/G NU (J+l ) 

0 =PL/GNU l 
XI =GL/GNU ( J) 

X2=QL/GNU(J+1 ) 

X =QL/GNU 1 
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C REFLECTION COEFFICIENTS 

AMOl r (2.5*C5(01))/(01*C3 (01)) 
GAM01=ATANCAM0t) 

AMX1 = <2.5*C5(X1))/(X1*C3(X1 )> 
GAMXl-ATAN(AMXl) 

AM02=(2.5*C5(02))/(02*C3 (02)) 
6AM02=ATAN(AM02) 

AMX2:(2.5*C5(X2))/(X2*C3 (X2)> 

GAMX2=A TA N(AMX2) 

DI FGM1 - 5 7.295 R* (GA MX I -GA MO 1 ) 

DIFGW2=5 7.295R*(GAMX2-GAM02> 

C0N= (EM*PERM*3.0ER) /(E* Efr ! ,5E3*5 7.295R) 

C CALCULATE ELECTRON CONCENTRATION 
SUB =<3L*C3 CX)-PL*C3 (0) 

GNU2=GNU1**2 

DELTA =DPHI ( J) -DI FGM2+DI FGM1 
ENUMrC ON* DELTA* GNU2 
EN(J)=ENUM/SUB 
EN(J)=EN(J)/1 .0ES 

50 WRITE(6,S0)HT,DPHI (J) ,GNU! ,EN(J) 

60 FORMAT(2X,F5.2,2X,F10.4,5X,E10.4,2X,F12.6) 

GO TO 10 
END 


